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We define a new observable that depends on finite redshift differences of the spin- 
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weighted angular moments of the two point function of the three dimensional cosmic 



Qh shear and on luminosity distance. It is shown that precise measurements of our 

^3 \ observable will be able to tightly constrain the expansion history of the universe. 
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The observational evidence that the expansion of the universe is being accelerated by 

dark energy (DE) is strengthening continuously 0, 0, U, [| . On the other hand, the theo- 

V/^ ! retical efforts to understand the nature of DE and the underlying reason for the accelerated 

OS ■ 

O expansion were less successful so far. It has become clear that in order to make further 

: 

theoretical progress better methods of determining the expansion history of the universe are 



needed. The kinematic distance probes in the homogeneous and isotropic universe, such 



as luminosity distance (di), angular distance etc., are limited because they rely on light 

& : n 

emitted from distant sources, and hence measure an integral over the expansion history [6J. 
It is possible to circumvent some of their limitations by adding prior assumptions on the 
evolution, or by combining several kinds of kinematic observations ODD , or by focusin, 
on better determined quantities A recent review of the subject can be found in 
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In addition to the kinematic probes one can use dynamical (fluctuations) probes to de- 
termine the expansion history of the universe. Cosmological perturbations provide, through 
their dependence on the homogeneous and isotropic background, an independent source of 



information 
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141 ]. One needs observables that can determine in a reliable way the time 
dependence of the perturbations' amplitude and not only their spectrum at a fixed redshift. 
Such observables can have a simple functional dependence on the expansion history and 
hence can provide much more information than a measurement of a single spectrum. How- 
ever, to provide significant constraints above and beyond those that kinematical probes put 
on the expansion history such observables have to be measured and evaluated precisely. 
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The three dimensional cosmic shear (3DCS) is a redshift dependent tensor that measures 
the shape distortion of distant galaxies as the light emitted from them propagates through 
the perturbed universe Q|. The two dimensional cosmic shear was recently measured 



in weak lensing (WL) experiments Il8l. Il9l| while planned programs are expected to 



have some three dimensional capabilities [20]. Several recent reviews of the theory and 



observations of WL provide an excellent description of the field |27|, |28(. Our results suggest 
that precise measurements of the 3DCS will significantly enhance our ability to constrain 
the expansion history of the universe. 

If the cosmological perturbations are small, the cosmic shear (CS) depends linearly on 
them. In this linear regime both observations and theoretical calculations can be done 
accurately. However, to take full advanta ge o f the linear regime one needs to filter the 



influence of the late time nonlinear regime 



29[. This requires the capability of accurately 



measuring finite differences of the 3DCS at different redshifts. 

The calculation of the 3DCS relies on solving the geodesic deviation equation in the 
presence of the cosmological perturbations to calculate the rate of divergence of light propa- 
gating on nearby geodesies. Obviously, this requires knowing the rate of divergence of light 
propagating on nearby geodesies in the absence of the perturbations, which gives the angular 
distance in the unperturbed universe. Hence accurate knowledge of the angular distance (or 
equivalently, di) is necessary for accurate knowledge of the 3DCS. 

In this paper we define a new observable J(z) derived from the 3DCS, and discuss how 
precise measurements of J can be used to tightly constrain the expansion history of the 
universe. The observable J(z) depends on finite redshift differences of the spin-weighted 
angular moments of the two point function of the 3DCS. The equations connecting the 
metric perturbations and the CS that allow us to make the quantitative estimates in this 
paper were developed in jscj. 

We consider the perturbed spatially flat Friedman- Robert son- Walker universe. We as- 
sume that the cosmic fluid is a perfect fluid and that the perturbations are adiabatic. In 
this case the line element of the perturbed universe is defined in terms of single space-time 
dependent perturbation $ 

ds 2 = a 2 { V ) [-(1 + $ (77, £)) d V 2 + (1 - $ (77, x)) (dw 2 + w 2 dQ 2 )} . (1) 

Here rj is conformal time and w is a radial coordinate. In [30] we have followed the standard 



derivations that are reviewed in |31[ to obtain the equations of motion for $ in terms of the 
total redshift dependent equation-of-state (EOS) parameter w tot (z) = p(z)/p(z). 

We briefly review the results of . When the universe is matter dominated the growing 
mode of $ factorizes <&+(z,x) = C(x)§t{z)- The equation for $t in terms of Wtot 

5 + 3w to t(^) 5 + Sw tot (z) 

can be solved as follows. First, define an integration factor 

1{z) = r 2(i+,) • (3) 

z 

Then reduce eq. (J2J) to a first order equation 

e-'Wft [e'«<M.)] = (4) 

whose solution is 

z 

The initial conditions are set at the initial redshift z- m . Asymptotically, for very large z the 
universe was matter dominated, hence w to t and $t is a constant whose value can be 
determine from eq.(J2J): $t — ► 1/5. 

We assume a power-law spectrum for $ + . In momentum space this implies that 
g 3 (C((fi), C*(g*2)) = ^4(<?/<?o) n-1 <5(<?i ~~ ^) where n is the spectral index and A is the spectral 
amplitude at go- Both were measured most recently by the WMAP and SDSS experiments 
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]. In particular, the spectral index is approximately n = .95 corresponding to an ap- 



proximately flat spectrum. The range over which the spectrum is flat (n is approximately 1) 



is limited because causal processes inside the horizon suppress the perturbations 



3aJ_p 
[30]) 



3J. (For 



further discussion see [30]). In the rest of the paper we will assume for simplicity a flat 



spectrum n = 1. This assumption can be relaxed 33]. 

The CS tensor 7 quantifies the distortion pattern of distant objects as light emitted from 
them propagates through the perturbed universe. It can be measured in WL observations. 
We will be interested in the angular moments of 7. Because 7 is a tensor it has to be 
expanded in tensor spherical harmonics 7 = J2im 2 a Zm iXimi and the two-point correlation 
function of 7 has a similar expression (7,7*) = 2C1 2^,-2- The relevant properties 



of the s- weighted spherical harmonics s Yi m can be found in [34]. The spin- weight operator 
5 operates on the s-weight spherical harmonic and gives an s + 1-weight spherical harmonic 

3 sYl,r~= [{I -s)(l + 8+ l)] 1 / 2 s+ lY ljm . 

In we have shown that 

-y(w,9,<p)= ldw ^ W ~ W > m$(w', 9,<p) (6) 

J W w' 





and that the angular coefficients 2 C\ are consequently 

2< ' "'" ,2 "" - 



U>1 

!Cl(wi , W2 ) = (teJMlfl + 2) / ^ (m '-" )(ro2 - M) j>| M , < «*. (7) 

J It Wi w 2 



We now define the new observable J (2;). As we have explained, to take full advantage of 
the linear regime the observable has to subtract the region of small u in eq. (f7j) where the 
integral formally diverges. However, the subtraction has to be of a kind that will allow a 
good measurement of the defined quantity. After considering different possibilities we have 
come up with the simple finite difference 

A W2 2 Ci{wi, w 2 ) = (^p^/ + 2 ) ^ ^Ciiw!, w 2 + Aw 2 ) - 2 Ci{w u w 2 ) 
Evaluating A W2 2 C[(wi,w 2 ) using eq.(J7|) we obtain 

A W2 2 Ci(w!,w 2 ) = — I — / du (wi -u) $%(u). (9) 

w x w 2 w 2 + Aw 2 J 


This suggests the definition 

w 

J(w) = I du{w-u) $ 2 T {u). (10) 





Since J(w) is I, w 2 and Au> 2 -independent it can be estimated by statistically averaging 
measurements for many different /'s and many different values of w 2 and Aw 2 , 

J(w) = Iw w 2 W2 ~^ w W2 A W2 2 Ci(w,w 2 ) \ , (11) 

\ ' l, W2, At02 

the only restriction is that w 2 > w. 
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Since w(z) = J = we may express J as a function of redshift 

W*55j(££-£SW (l2) 

v / 

We would like to show that accurate measurements of J can significantly enhance our 
ability to constrain the expansion history of the universe. We will use eq. flllj) to determine 
the requirements for an accurate measurement of J. First one has to measure accurately 
the angular moments 2Ci(w\, W2) for many Vs. Since the measurements of the CS are noisy 
because of intrinsic ellipticity variations and observational limitations one uses statistical 
methods to reduce the errors. The more galaxies are measured, the better the shear mea- 
surement. The best that one can hope to do is to measure the angular moments such that the 
errors are dominated by statistics. Hence, the ultimate limitation is set by cosmic variance. 
By averaging over many angular moments from l\ to 1% 1% ^> 1) one can reach a relative 
accuracy of ~ \j \Jl\ — l\ . The largest I for which the 2C1S can be measured accurately 
and the averaging can be done reliably will determine the ultimate accuracy. Thus, to get 
below 1% accuracy one would need to make a cosmic variance limited measurement of about 
a 100 angular moments in the range I ~ 100. 

Precise measurements of the angular moments are not sufficient. In addition w(z) and 
Aw(z) have to be measured precisely. This means that the redshift at which the angular 
moments are measured and for this redshift have to be known precisely. If one wishes 
to obtain a 1% measurement of J then redshift measurements and di measurement at or 
below 1% have to be made. While a 1% redshift measurement is significantly below the 
expected accuracy of photometric redshift measurement it can, of course, be measured to 
much better accuracy using spectroscopic methods. A 1% measurement of di is a plausible 
goal for future supernova observations, however, accuracies of much below 1% are not likely. 
Hence the final limiting factor on the accuracy of measuring J is likely to be set by the 
accuracy of di measurements. To summarize, measuring J with one percent precision (or 
perhaps somewhat better) seems a plausible optimistic goal, though a great challenge for 
future observations. 

We wish to estimate the sensitivity 5J/5w tot . It is very instructive to obtain first the 
analytical the dependence of the sensitivity on different parameters. We start by calculating 
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S^T/^w tot . From eq.Q we find after some algebra that 

^ ( *t{z) jdz + 5* T (zU - e™*M«) = (13) 



e 



dz- 



l + w tot (z)) I dx y^*? X \ + Sw tot(z) 
2(1 + x) 



2(1 + 5) 

z 

Equation (|T3*j) determines in terms of 6w to t- 

{€ Zi < Z < Z 2 
. In 
Otherwise 

this case we may choose the upper limit to be Z{ = z 2 . To calculate 5§/5w tot it is useful 
to consider the limit z 1 2 ~ z zi <C 1 and to assume that $(zi) is fixed so its variation vanishes. 
Additionally, $t ~ 1/5 in our normalization. In this case we obtain s -§^ ~ e g^y. The 
relative change increases linearly, the maximal relative change being at zi, ig^r^y — e (i+zl) • 
The relative change is linear in the absolute change e in w to t and linear in the width Z\ — z 2 
where the change occurs. Earlier changes in w tot produce smaller changes in $t- We now 
turn to the range z < z\. Recall that for z < z\ w tot is unchanged, so it is better to choose 
Zi = Z\ in eq. (|13|) . In this case 51 = and the only source of change is now the initial variation 
of $ r (zi). Additionally, I(z x ) = 0. It follows from eq.flHU) that 5$> T (z) = e- J ^5$ T (zi). To 
get an idea about the numerical factors consider the case that Wtot is approximately constant, 

5+3w tot 5+3m tot 

thene-^(^) 2 , and ~ e 2 . Since 5/2 > ^ > 3/2 

for z < Zi, 5&t/$t slowly decreases from its maximal value at Z\. Our estimates are in 
good agreement with the numerical results. 



We can now estimate 5 J, 5J(z) = ] dy ^ ( %g - ^ ] § 2 T {y) 



5H(y) Sd L (z) 



H(y) 1 d L {z) 



d L {y) " rz *t(s/) 



Again, this fits well with the numerical estimates. We should emphasize an 

important feature. If 8w tot vanishes for values less than some redshift Z* then 8di and SH 
vanish because they depend on an integral from z = to z*. In this case only <5$t is non- 
vanishing. For z > z* the situation reverses. Now 5$^ vanishes, and the only contribution 
comes from the other terms. In Fig. ^ we show the results of the numerical evaluation of the 
induced relative differences in J due to modifications of u>tot of the form discussed above. 

If the content of the universe is parametrized in terms of a cold matter component with 
relative density today of Q m and a DE component with a constant EOS parameter wq 
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FIG. 1: The relative changes in J(z) induced by modifications in u^ot- In the top panel we show 
the EOS of the £l m = 0.3, ACDM fiducial model and the EOS of three modified models. In the 
middle panel the relative deviations of the three modified EOS from the EOS of the fiducial model 
are shown, and in the lower panel the corresponding relative deviations of J. 

and relative density Qq = 1 — Q m , then the constraints on the function w to t{z) become 
constraints on the parameters Q m and wq. We discuss this example mainly to compare with 
the sensitivity of di measurements as is often done in the literature. 

In Figure El we show 1% and 2% distinguishability contours of J (filled) and d^ (empty) 
in the range z < 1.5. The fiducial model is a Q m = 0.3 ACDM model. Two models are 
deemed distinguishable at 1% if their J's (g?l's) deviate by more than 1% in the whole of the 
redshift range. Also shown are relative deviations of J(z) from the fiducial model for four 
different constant EOS models. The d^ contours in Fig. El represent much better accuracy 
than current measurements . They are approximately equal to the final precision goals 
of the Supernovae Legacy Survey la ~ 1%, 2a ~ 2% 35j. As can be seen from the figure 
a measurement of J with the same precision of a measurement of di gives about an order 
of magnitude improvement in determining the cosmological parameters. As emphasized 
previously, the figure should not be interpreted as a measurement of wq or of Q m , rather 
a measurement of the total EOS w tot . Despite the fact that the observable J is related to 
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FIG. 2: The upper figure shows 1% and 2% distinguishability contours of J (filled) and (empty) 
in the range z < 1.5. The fiducial model is a Q m = 0.3 ACDM model. In the lower figure we show 
relative deviations of J(z) from the fiducial model for four different constant EOS models labelled 
by (n m ,w Q ): (from top to bottom) (.32, -1.05), (.31, -1.025), (.29, -.975), (.28, -.95). 

perturbations in the matter only, and despite the fact that the J contours in Fig. El appear 
to be (1l contours with a limited range of fl m , Q m is not really measured separately here. 

We have used a simplified model to obtain our sensitivity estimates to which many im- 
provements and refinements can be made. We have found that obtaining a limiting accuracy 
as small as 10% over a range of redshift as small as 0.1 for local variations in io tot seems 
plausible. The implementation of our method will require extraordinary observational efforts 
because it requires very accurate measurements of the shear, very accurate determination of 
redshifts of the sources, and in addition very accurate measurements of as a function of 
redshift. To achieve such ambitious observational goals would require combining the results 
of several accurate experiments. 

We believe that determining the nature of dark energy is still an outstanding problem. 
An independent determination of Q m (if it is possible to do so) will turn the constraints 
on u>tot into a constraints on the EOS of the DE. However, there remain degeneracies for 
groups of expansion histories that oscillate about each other [8(. This type of degeneracy 
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can perhaps be overcome by combining di and other observables with measurements of J. 
In any case, we do expect some remaining degeneracy that will have to be resolved by better 
theoretical understanding of the nature of dark energy. 

This research was supported in part by the National Science Foundation under Grant 
No. PHY99-07949. 
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